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1 CODED LABEL INFORMATION EXTRACTION METHOD 



2 Introduction 

In certain types of human and machine-readable information-bearing labels a set of 
elements is used to represent the information contained in the label. This 
representation may be by varying the characteristics of the elements comprising the 
label, and also by the position in which the label elements are placed. Reading 
apparatus senses the characteristics and placement of the elements in order to 
decode the information contained within the label. The elements in the label are 
sequentially scanned, to discover the presence of die constituent elements and to 
measure their characteristics and position and thence to decode the information 
contained by the label. 

This appUcation describes label decoding methods by which information represented 
on such a label can be decoded. Typical label embodiments are where the label is 
manufactured from a plurality of magnetically active elements supported on a 
substrate. Information is coded by controlling the relative positions. Independent 
control of the physical properties of the elements, such as shape may also be used. 

In alternative embodiments, information may be coded by controlling die relative 
angles of orientation of the easy axis of the elements, in some cases together with 
independent control of element shape, or otiier physical property. 

3 Description OF THE Invention 



3.1 Scanning Methiod 

In certain types of label decoding system a label consisting of a set of elements is 
scanned, that is the elements are activated in turn by appljdng a signal to which the 
scanned element responds, and to which the reading apparatus is sensitive. Either 
or both of the means of generating the scanning signal, and the means of generating 
the received signal in the reading apparatus select a single elements in turn so that 
the signal from a particular element may be resolved from that of die other elements. 

The means of scanningmay be automatic by electronically scanningthe interrogating 
signal and the receive sensitivity along die label or by an automatic mechanism 
causing the label elements to be interrogated in turn. Alternatively the scanning may 
be manually operated where the label is moved past the reading apparatus by hand, 
or alternatively the reading apparatus or part thereof may be moved past the label 
by hand. 

3.2 Signal received from reader 
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Fig la shows a plan view of a label where the elements are made out of high 
permeability, low coerdvity magnetic thin film material such as Atalante 
manufactured by IST(Belgium). In this case die elements are lmmxl-5mm. The easy 
axis of the material is arranged to be parallel widi die longitudinal axis of the label. 

The reading antenna is passed by the label subjecting each element in turn to a 
region of null magnetic field formed by a permanent magnet arrangement. A 
superimposed high frequency field, in this case 6.348 kHz causes the material 
element to switch when it is present in the magnetic null. A receive coil is used to 
detect the harmonic signals resultuig from the element switching, and die received 
harmonic signal is mixed to baseband and filtered prior to digitisation. Typical 
signals received from such a label, detecting the second harmonic of the element 
switching signals, when the read head was passed in proxunity to the tag are shown 
in Fig lb. The signal from a single element is shown in Fig Ic This shows two half 
peaks of opposite sign with a zero crossing. It can be seen that the label signal 
approximates a superposition of single element responses at each position m which 
there is an element 

The signal is sampled witii a fixed sample rate, and converted to a digital 
representation in an analogue to digital converter. The sample rate must be fast 
enough to capture the signal information in accordance with the Nyquist criterion, 
preferably with anti-aliasing filtering prior to the sampling operation. 

The sample nmnber of the centre of each element response depends upon the time 
at which the element passed the read head. 

3.3 Velocity estimates from signal 

In addition, the width of the element response depends upon die speed at which die 
element passed the read head. Information about the dynamics of die relative motion 
of the tag and die read head can be inferred from measuring the width of die signal. 
This can be achieved by measuring the width of the positive going and negative gomg 
peak associated with the element. It can also be achieved by measuring the 
separation between the positive and negative going peaks. This speed estimate is 
essential to estimate the dynamics of motion and to decode the information in the 
label reliably without compromising the available code space by adding regidar 
features in the label. 

In alternative implementations of die reading apparatus, the scan n ing motion of the 
magnetic null relative to the label may be achieved by electronically scanning using 
a set of excitation coils that are driven in a phased manner so as to cause the 
magnetic null to propagate in space. The null may be scanned at sufficient speed that 
the switching of the label elements may be detected direcdy from the receive coUs 
without needing to superimpose any separate high frequency excitation. In this case 
die material element response coupled into the receiver is a single peak. Fig Id 
shows 3 adjacent peaks of the fundamental signal from 3 magnetic elements 
recorded from areading system employing electronic scanning. In this case the width 
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of the response from each magnetic element can be used to provide ^^^^^^^f 
Ae velocity of the scanning null as it passed over the element bemg scanned. This 
S^ormaSJ^ can be used to^estimate the detaUed dynamics and 
SdXo from misaUgmnentof the label from the nomuial read axis of the read head. 

3.4 Amplitude Coded Signak 

Information can be encoded into the label by varying the height of ^f^^^^ 
Semite, thatis thesizeofthedementnormaltothelonptudindaMsof thelabel, 

or Sr^Sdth, that is the size of the element along the longitudmal a«s. 

By mamtahung the element width constant, variations in ^^^^ height 
pLdoSantlyLuseachangeinampUtudeofthereceivedsignalasiUu^^^ 
2rfor the mechanicaUy scanned second harmonic system, and m Fig 2b for an 
^e^oSX^ed^rstem detecting tiie fimdamental signal received from the 

^ label. 

A e«t «f l,n«.ar test labels with 3 elements was manufactured with the element 
posfiol^S^ i J^^^^ 

SSdS^ wie imii and the height of the centre element varied bet«reea 
^Tv^^ctx label, m Fig 2a. the peak position, tiie half peak width, the half peak 
s^ation^Setiph^dedifferL^ 
25 Sve v^^tion of thLe parameters of the central test dement -ompared^^^- 

outo elements, which were mamtained at constant hei^t. were plotted The 
rSurlrcSS^lelementampUtudeCvampr) relativ^^^^^ 
a dear, near Unear relationship to central element height The 
dS^St^elative position (vposr). widfli (vwidr) ^^/^P^^^^^^f^^'^P^^^f 
70 do not vary significantiy witii central element height, over and above tiie effects ot 

'° ^^SS^rs and unperfections m tiie manufacture of tiie test labels. 

Similarly Fig 2b shows tiie relative variation of measured centi^ element ampUtude. 

IZiZn aSi widtii witii element height when tiie test labek w«-e measured m an 
• r<So^Sly ^Lmmg reading system and shows a shnilar near Unear f elationship 

^ to S mechmiically s<Lned readmg metiiod receivhig second harmomc signals. 

Element height can be used to code additional mformation witii a detector based on 
W^deTtiie receive processing. In flie simplest case, two hf^^ts are used to 

40 ^d^^JtrelmentsinLlabel.TMsismustratedinFig3^^^^ 

witii elements of one of two heights. The correspondmg '^^f^^^^"^'^^^^ 
Fig 3b. The difference m received ampUtude from elements of different h«ghte can 
Kenned. A typical design goal is to maximise tiie ^^f^^^Jl^^^^^ll^, 
read A compromise is needed to maximise readmg range for a 2;height lab^, at me 

4. s^e tiSrS maxunismg die abihty to discriminate tiie heights of mdmdual 

eJle^ to order to discriminate die height of label elements rehably, certam of 
tiS^toSite may be constrained to be at a known height to act as a ref«-^ce hi 
^e e^lTiUus Jated m 3a, die 2 end elements of die label ^'^^^^^f^ 
^^^S^hei^t to act as a reference. A design rule may also constrmn a certain 
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number of the inner elements to also be at maximum height The reference elements 
may be at any a priori fixed position, or at a position determined by an a priori rule. 

The half peak width, or half peak separation can be estimated and used to form an 
estimate of the velocity of the scanmng, at the time when the element was scanned 
in the case of a reader detecting the second harmonic of the label signal' 
Alternatively, m the case of a reader detectmg the fundamental of the label signal* 
an estimate of the peak widtfi can be used to form an estunate of the velocity of the 
scannm& at the time when the element was scanned. These estunates are lypicaUy 
K/vf_est where K is some constant, and w_est is the estimate of width or half peak 
separation, as appropriate. The value of K typically depends on parameters of the 
reader, the element shape and the position of the element relative to the read 
antenna. Smce K varies with element properties, to facilitate estunation of the 
scannmg dynamics, the process to estimate K may be restricted to elements of the 
same type based on the amphtude of the element signals. 



3.5 Peak finding algoiidun 

The first step m processing a received signal such as shown m Fig lb or Fig 3b is to 
J^f parameters of the peaks contained withm the signal. This is achieved 
with the following steps: 

1 The signal is delimited, that is the continuous received signal is broken up mto 
a contiguous set of samples contammg the received signal from one label, or a 
set of peaks that is likely to contam enough mformation to decode the label m 
the case of non-registered repeating codes. ' 

2 The signals from this delimited period are sorted mto order of ascendmg 
mapitude. An indexed sort function, for example using tiie Quicksort metiiod, 
such as mdexxO (Numerical Redpies m C, second edition, ISBN 0 521 43108 5) 
j^^P»'0Priat« so tiiat tiie original location of the point of any rank can be quickly 

3 The positive half peaks are found by workmg down from die largest signal m 
magnitude m die sorted data, testing die signal points as follows. If the test signal 
pomt is a peak, that is its adjacent signals are lower m magnitude it is estimated 
as a positive half peak, and its amphtude is hiterpolated by interpolation of the 
peak and its neighbouring points. 

4 From each positive half peak found, the signal is scanned forward and backwards 
to find tiie points at which it faUs below a defined magnitude relative to tiie peak, 
this is interpolated to a fraction of a sample period A typical magnitude value is 
0.8 of the peak value found. The half peak widtii is estimated as die difference in 
position between diese points. The half peak position is estimated as tiie average 
of diese two points. The half peak is rejected if it overlaps witii any otiier pre- 
ensting peak m die Ust, tiiat is its widtii about its position interferes tiie widdi 
about any odier peak position in die Ust. If it is not rejected, die half peak is 
appended to a list of half peaks. 

i The peakfindingprocess stops when die magnitude of die next highest pomt falls 
below an absolute du-eshold, or when it falls below a defined fraction of die 
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5 largest magnitude signal, for example 0.3 of the largest magnitude signal, or if 

a limit for the number of peaks is reached. 

6 The negative half peaks are found by working up in magnitude from the lowest 
signal in magnitude in the sorted data, testing signal points of increasing 
magnitude. If the test signal point is a negative peak, that is its adjacent signals 

1 0 are hi^er in magnitude it is estimated as a negative half peak, and its amplitude 

is interpolated by interpolation of the peak and its neighbouring points. 

7 For each negative half peak, the signal is scanned forward and backwards to find 
the points at which it increases above a defined magnitude relative to die 
negative peak, this is interpolated to a fraction of a sample period. The half peak 

15 width is estimated as the difference m position between Aese points. The half 

peak position is estimated as the average of these two points. The half pe^ is 
rejected if it overlaps with any other pre-existing peak in the list; Aat is its width 
about its position interferes the width about any other peak position in the list 

^ If it is not rejected by this test it is then appended to a list of half peaks. 

20 8 The negative peak finding process stops when die magnitude ofthe next highest 

point falls above an absolute threshold, or when it falls above a defined fraction 
of the largest magnitude signal, or if a limit for the number of peaks is reached. 
9 The resultant list of positive and negative half peaks is then sorted into range 
order. This process will interieave the consecutive half peaks. 

25 10 Merging range consecutive pairs of half peaks in the range sorted half peak list 

forms a list of double peaks. The amplitude of the merged double peak is the 
di£feren<% in amplitude of two half peaks being merged. The position of the 
merged double peak is the average of the two half peak positions. The width of 
the merged double peak is the average of two half-peak widths. The separation 

30 of the merged double peakis the difference between the two half peak positions. 

The results of diis half peak finding process are illustrated in Fig 4a which shows the 
half peaks found when the signal shown m Fig 3b was processed, and shows the 
output of step 9 above, these peaks are shown Widi a cross at the located position 

1^ and estimated ampUtude. Fig 4b shows the magnitude of the double peaks in 

amplitude units. In this case tide amplitudes are negative since die first half peak is 
a negative going peak. The difference in peak amplitude between the full hei^t and 
reduced height elements in the label can be dearly seen. Fig 4c shows the double 
peak separations versus double peak position showing the variation across the 

40 scanning of the label. 

In the case of a fundamental signal where there is just one peak per element; such 
as that shown in Fig Id, a simpler process is may be used, missmg steps 6 to 10 
above since there is no need to find and merge double peaks. 

45 10.1 Extraction of velocity function 

There are N label elements positions whose centres are located at X, .. along the 
longitudinal axis of the label, and whose position we wish to estimate. It is assumed 
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that Xj is at position 0.0 and that the final element, X^^, is at the known position of 
corresponding to the length of the label L. 

The peak list as described above contains estimates of the sample number, and hence 
times Tj at which die reader detected a double peak. For each element i, the time T,. 
is simply the position in samples where the peak was measured multiplied by the 
sample period, which is the reciprocal of the sampling rate. 

By using the estimated double peak width or double peak separation for each of the 
N peaks denoted Wj a velocity estimate for the i* peak is VE,=1/W, at time T, . 

The order M velocity function of time is modelled by a polynomial function. M is 
chosen to allow sufficient degrees of freedom to moddl the dynamics of label motion 
sufficientiy accurately. It can range from 1 (constant acceleration model) upwards 

VF(t,M)=ao+ait+a2t^+a3t*+....aMt*' 

where {a} are a set of polynomial coefficients to be determined* 

The squared error for velocity estimate VEj at time T, is Ei*=(V]^-VF(Ti,M))* 

The total sum of the squared errors is ES(M)=Ei^+ + E,^^ which is 

minimised by varying the M+ 1 parameters a^o.^ a^ This is achieved by the method 

of least squares to find the optimum parameters in tiie least squares sense. 

Define the matrix X to contain the T| values raised to a given power. For example 
coliunn 0 row i of X contains and coliunn M row 4of X contains T4M . Define tiixe 
vector V to contain the velocity estimates VEi . For example tiie element in row 4 of 
V contains VE4 

Some linear combination of tiie columns in X will be closest to the vector V in the 
least-squares sense. The coefficients of this linear combination will be the 
coefficients of the best-fitting polynomial. If A is the vector of length M + 1 that holds 
the coefficients, the problem is to find the value for A that gives the best "near- 
solution" to the system of equations XA= V. If the columns of X are independent, the 
best choice for A is given by: 

A=(X^X)-'(X^V) 

Where represents matrix transpose and •* " represents matrix inversion. The 
coefficient a^ is then the element from row i of A. This may be effidentiy calculated 
by using methods such as LU decomposition, for examples the functions ludcmpQ 
and lubksbO (Numerical Recipes in C, second edition, ISBN 0 521 43108 5) 

This process is illustrate in Fig 4d where the velocity estimates that were 
computed from the double peak separation data of Fig 4c are shown, those points 
marked ''full data** contain the velocity estimates from all the peaks. In this example 
case, only the estimates from peaks with high amplitude, see Fig 4b, are selected for 
velocity fitting. These selected points are marked ''Data points used in fit". If the label 
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known a priori to have elements of the same size then all elements would Imve be«i 
used. A least squares fit of order 2 was used and the continuous hne shows the 
estunated velocity function as a function of sample number. A least squar^fit of 
order 1, 2 or 3 is typicaUy sufficient to model the velocity function with suffiaent 
accuracy. 

Error metrics are computed to allow the decoding algorithms to assess the quality of 
the read. A useful normalised velocity error metric, denoted evel, to quantify how 
wett the velocity function fitted the data observed is defined. It is the sum of squares 
of die normalised error definied as : 

evel=sum( E^*/VE^*) fori=l..N 

The vector A could also be estimated by an iterative minimisation procedure such as 
Powell's method. 

The above procedure can be repeated for varying model order M on the same data 
in order to find the best approximation. Preferably this is the solution wiA the 
smallestnumber of degrees of freedom, that fits the data read by the reader with the 
lowest error metrics. 

10.2 Analytical integration and normalisation of tag length 

The velocity fit function is VF(t,M)=ao+a,t+aat?+a,t?+....+a„t" This may be 
analytically mtegrated to find the positional fit function 

P(t,M) =Q+aot+(l/2)ait2+(l/3)a2t'+....+(l/(M+l))a„t"** 
where Q is an arbitrary constant of integration, set to zero. 

We know that the position of label element 1 is at 0.0. We also know Aat theposition 
of element N, the last element scanned is L, die lengdi of the label. The final 
estimates of element position are: 

XE, = L . (PCX,, M)-P(Xi,M)) / (PCXn,M)-PCXi,M)) 

10.3 Binning of data 

There are N positional estimates XE^ to XE^ from which N-1 gaps may be estimated. 

Typically the gaps have the minimum gap subtracted, and are divided by Ae 
minimum gap, to givearealnumber representing thenumberof gap uioremCTite.'nie 

nearestmteger is then taken to representthe gap for subsequent decodmg. This can 

be e3q>ressed as: 

The N-1 normaUsed gaps are calculated by r _ i m i 

Rj=(XE,+,-XE,-Minimum_Gap)/GapJncrement forj-l...N-l 
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The normalised gaps are tben converted to their nearest integer representaton 
Gj= int(R|+0*5) for j=l...N-l, whereint(x) means the truncated integer portion 
of X. 



One useful error metric, denoted emax, is defined as the maximum magnitude of the 
difference between Gj and R, Mother useful error metric^ denoted eav, is the average 
magnitude of the difference between Gj and R|. 

Fig 4e shows the gap sequence data R| computed from the data in Fig 4a-d, and 
shows the dose correspondence or "snapping^ to an integer grid, giving a high 
confidence that the gaps have been correcdy meastu-ed. The gap sequence in this case 
is 3,2,2,1,7,0,0,2,1 • 



10.4 Enumeration algorithms 

Depending on the coding scheme used, the gap sequence is enumerated to a unique 
decoded value. Various codes exist that have different properties for example those 
codes allowing read direction to be determined, or those that allow non-registered 
operation, that is allowing the gap sequence to start at any position, for example a 
portion of a repeating gap sequence. The code may be determined a priori, or 
diagnosed by examining the read signals for the ntunber of material elements 
present, or by estimating the sum of the gap ntunbers. It is preferable to determine 
this beforehand in order to minimise the chance of misinterpreting a given set of read 
data. 

For certain codes, the number of elements in the label may be one of a range of 
possible values, with fixed label length in order to increase the munber of possible 
codes. 



The direction finding codes are preferably used where amplitude coding is used to 
allow the amplitude information to be associated with the correct elemenL 

10.5 Amplitude threshold estimation of elements 

Where it is desired to detect the size of the received element signal, for example 
when the elements may be of two or more sizes, reference elements may be placed 
in the set of elements used to make the label. In one embodiment, the end elements 
are constrained to be die largest element size to act as a reference. Other 
embodiments may place reference elements that have some unique characteristic, or 
may be identified by an a priori rule. 

The amplitudes are measured and used to compute normalised amplitudes. This is 
because the reference elements maybe significandy different in magnitude. This may 
be due to the scanning process varying the distance or angle of the read head from 
the label during the read, or other factors. 
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The estimates of element position XE^ may be used for linear interpolation of the 
reference amplitude along the label. The reference amplitude predicted at the 
estimated position of the element divides the measured element amplitudes to 
normalise the amplitudes. The normalised amplitudes should be nominally 1.0 for 
an element the same size as the reference element. 

Threshold levels are set to provide discrimination between the possible element 
amplitudes and to classify the amplitude of the element into one of an alphabet of 
values. For example if diere are two values then die amplitude is either high or low. 
For example a half height element may have an amplitude, on average, half the 
reference amplitude ie 0.5. In this case, the threshold level should be set to 
approximately 0.75 in order to determine if the element was full heis^t or half 
heiglit. 

Preferably a code is used which enables the direction of reading of tiie label to be 
determined by the sequence of gaps that resulted. In this case the elements can be 
uniquely identified and it is known which element of the label to associate with the 
decoded amplitude information. 

In the case of the label data shown in Fig 4b, denoting the amplitudes of the signals 
from the elements as '*H" for high amplitude and "L** for low amplitude, the elements 
are categorised as H,L,H,L,H,L,H,L,L,H. A larger alphabet of elements could be used 
by using a label with more possible shapes, for example. The categorisations are then 
converted to decoded data. 

A typical scheme allows the inner elements only to contain data, ignoring the end 
elements, which are references, and maps H to a logical 0 and L to a logical 1 so that 
the binary data word from the example data would be 10101011. In tiiis scheme up 
to 256 amplitude codes can be represented from 8 inner elements. Preferably the 
codes used are constrained to ensure that a certain number of elements in the label, 
for example 4 are at full height to facilitate velocity estimation. If this constraint is 
used, the number of amplitude codes is reduced by 37 to 219 in this case. 

An ampUtude error metric eamp is defined which is the normalised amplitude 
difference between the largest half height element diagnosed to the smallest full 
height element diagnosed. If there is no half-height element diagnosed, eamp is 
defined as the smallest fidl height element amplitude. In a label containing different 
sized elements, a large value for this metric suggests a dear discrimination in 
amplitude between elements of different types. Smaller values indicate increasing 
risk that an element may be categorised as the wrong amplitude. 



10.6 Error Metrics 

The error metrics defined: evel, emax, eay and eamp can be used to determine the 
read quality, evel, emax and eav are preferably small and eamp high. Independent 
thresholds for these quantities are independently established. These thresholds are 
based on measured values for many labels in many read operations over the range 
of reading ranges and scanning dynamics expected in an application. This is done to 
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minimise the probability of a wrong code at the same time as maintaining a hi^ 
probability of a correct read. 

During the read decoding process, if the computed error metrics are hi^^er than the 
threshold value in the case of evel, emax and eav then the read is rejected. . If the 
error metric is lower than the threshold in the case of eamp then the read is rejected. 
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Fig lb: Typical received second harmonic signal from label 
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Fig Ic: Typical received second harmonic signal from element of label 
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Fig 2a: Variation of amplitude with element height^ second harmonic signal 
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Fig 2b: Variation of received amplitude with element height - fundamental 
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Fig 3a: linear label with height coding 
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Fig 3b: Receive signal from linear label with height coding 



G:\AMENDED\MJA-FN29 • coded label information extraction mediod l.wpd 



700 




p,peakpos 
Sunple muiiber 



Fig 4a: Detected half peak positions 
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Fig 4b: Double peak amplitudes versus time 

Element amplitudes categorised : H,L,H,L,H,L,H,L,L,H 



^746408^ 3 



^ dpcaksep 



§ o 



J.90719^ 




0 

^6.186223^ 



SO 



100 
dpeakpos 
Double peak position 

Fig 4c: Double peak separation versus time 
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Fig 4d: Velocity lit to selected points versus time 
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Fig 4e.^ Resultant fit of Gap nmnber to integer grid 

Nearest integers to gap number is 3,2,2, 7,0^0^2^1 
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